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Abstract 

We assert that the physics underlying the extraordinary light transmission (reflec- 
tion) in nanostructured materials can be understood from rather general principles 
based on the formal scattering theory developed in quantum mechanics. The Maxwell 
equations in passive (dispersive and absorptive) linear media are written in the form of 
the Schrodinger equation to which the quantum mechanical resonant scattering theory 
(the Lippmann-Schwinger formalism) is applied. It is demonstrated that the existence 
of long-lived quasistationary eigenstates of the effective Hamiltonian for the Maxwell 
theory naturally explains the extraordinary transmission properties observed in various 
nanostructured materials. Such states correspond to quasistationary electromagnetic 
modes trapped in the scattering structure. Our general approach is also illustrated with 
an example of the zero-order transmission of the TE-polarized light through a metal- 
dielectric grating structure. Here a direct on-the-grid solution of the time-dependent 
Maxwell equations demonstrates the significance of resonances (or trapped modes) for 
extraordinary light transmission. 



PACS: 42.79.Dj, 41.20.Jb, 42.25.Fx, 03.65.Nk 
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1 Introduction 



Supported by technological progress, studies of the interaction of electromagnetic radiation 
with nanostructured materials have become an area of intense research driven by potential 
applications in optics and photonics [l]-[4]. In particular, it has been found that metal [3]- 
[16] and dielectric [17]- [24] grating structures posses extraordinary transmission (reflection) 
properties within narrow intervals of wavelengths close to the grating period. While for 
dielectric gratings a common point of view on this phenomenon, as occurring due to the 
existence of trapped modes or guided wave resonances [18] -[24], seems to be established, there 
is still an ongoing discussion on the origin of a nearly 100% light transmission within narrow 
wavelength range (s) observed in slit and hole arrays made in metal films. Although similar 
results are obtained with different theoretical approaches, in which the Maxwell equations are 
numerically solved, an explanation of the underlying physics, as due to excitations of coupled 
surface plasmons and/or cavity modes [5, 6], competes with the dynamical diffraction theory 
point of view [9, 10]. Remarkably, all available studies of metal grating structures have been 
carried out with the TM-polarized light (the magnetic field is parallel to the slits) where 
surface plasmons can indeed be excited. 

Transmission and reflection properties of grating structures are typically studied by sta- 
tionary methods in the frequency domain. Nevertheless, the dynamics of light scattering 
on gratings can partly be guessed from stationary results. Indeed, consider a femtosecond 
(broad band) pulse impinging on a grating structure such that the spectral range of the pulse 
is much larger than the grating transmission window(s). From the uncertainty principle it 
follows that, in order for transmission of light to occur only within a narrow frequency range, 
the radiation should last much longer than the duration of the initial pulse. This implies that 
the corresponding electromagnetic modes have to be trapped by the nanostructured mate- 
rial during a sufficiently long time. Such modes are known in scattering theory as scattering 
resonances [25]. 

Here we propose a general point of view on the physics underlying the extraordinary 
light transmission in nanostructured materials. By reformulating the Maxwell equations 
in the form of the Schrodinger equation and by using quantum scattering theory, we show 
that this phenomenon is a direct consequence of the existence of trapped electromagnetic 
modes possessing large lifetimes. Based on this general principle, various structures can be 
custom-designed that would transmit or reflect light within a designated narrow wavelength 
range(s). As an example to illustrate our approach, transmission properties of a metal- 
dielectric grating structure are studied for TE-polarized light (the electric field is parallel 
to the slits) impinging normally on the grating. A direct on-the-grid solution of the time- 
dependent Maxwell equations demonstrates the significance of resonances (or trapped modes) 
for the enhanced light transmission. Interestingly enough, such a system has never been 
studied before, probably because of the absence of the coupling of electromagnetic modes to 
plasmons, often thought to be the main mechanism of the extraordinary light transmission. 
It should be understood that our choice of the geometry does not imply an attempt to favor 
either the plasmon or cavity-mode point of view in, generally, strongly coupled systems. 
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2 Metal— Dielectric Grating 



We begin with an example of a grating structure sketched in Fig. 1 to illustrate our basic 
idea. The grating structure has translational symmetry along one of the Euclidean axes, 
chosen to be the y axis. The structure is periodic along the x axis with period D g , and the z 
direction is transverse to the structure. For the sake of comparison with previous works [5], 
the parameters are chosen to be: the grating period D g = 1.75 fim, the thickness (along the 
z axis) h = 0.8 -j- 1.4 fim, and the opening width a = 0.3 -j- 0.35 fim. In the case of the TM 
radiation impinging on the grating with no dielectric fillings, the extraordinary transmission 
properties have already been reported and analyzed [5, 7, 10, 11, 13]. The question arises 
whether enhanced transmission can be obtained for the TE radiation, or according to our 
remarks in the previous section, whether or not trapped modes exist and can be excited. 

To obtain a rough estimate of the wavelengths of possible trapped (or quasi- stationary) 
electromagnetic modes in the system, consider first the case of a perfect metal grating. For 
a moment, we also neglect the effects due to a finite thickness of the grated metal slab in the 
z-direction. In other words, we neglect the coupling between trapped modes and radiation 
modes outside the grating. Then the x component of the wave vector is quantized as a 
consequence of the zero boundary conditions at the metal-dielectric interface. A quantization 
of the z component of the wave vector can be understood as Fabry-Perrot modes in a 
dielectric slab inside the metallic waveguide, i.e., the modes for which the dielectric slab 
is 100% transparent. Admissible (quantized) values of the wave vector inside the grating 
define wavelengths \ nm of trapped electromagnetic modes that can be excited by the incident 
radiation 



where e is the dielectric constant, and n,m — 1, 2, .... Note that, since a < h for the grating 
geometry, the first term in the r.h.s. of Eq. (1) is the most relevant for the threshold to 
excite quasistationary states (standing waves). In what follows we are interested in the zero 
diffraction mode for wavelengths A > D g so that the reflected and transmitted beams prop- 
agate along the z-axis. In this case the threshold for trapped modes to exist is determined 
by the lowest mode, n — m — 1, leading to the condition 



Taking into account the grating parameters, condition (2) can be satisfied only if e > 9.7 for 
the range of a and h specified above. Thus, if there is no dielectric filling (e = 1), the grating 
structure would totally reflect the TE radiation with A > D g , which is indeed supported 
by our numerical simulations. This is in contrast to the TM radiation case where 100% 
transmission can be reached for the very same grating with no dielectric fillings and made 
of the perfect conductor [5] . 

In our numerical study, we consider the grating with openings filled with a non- dispersive 
dielectric with e = 11.9 (this corresponds to Si within the wavelength range under consider- 




(1) 
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ation). The metal is described by the Drude model: 

u 2 

e M {u) = 1 - p . . (3) 

For the sake of comparison with previous works the plasma frequency and the attenuation are 
taken as: uo p = 9eV, and 7 = O.leV [5]. Transmission and reflection properties of the grating 
structure are calculated by means of the wave packet propagation method. The method is 
based on the representation of the Maxwell equations in the form of the Schrodinger equation 
for which the initial value problem is numerically solved by a time-stepping algorithm. A 
detailed description can be found elsewhere [23, 24]. 

For the specific case considered here, the Hamiltonian formalism is as follows. Let D = 
E + P where D, E, and P are the electric induction, the electric field, and the medium 
polarization vector, respectively. In the Drude model the medium polarization vector satisfies 
the second order differential equation 

P + 7 P = ^ 2 E, (4) 

where derivatives with respect to time t are denoted by overdots. Equation 4 must be solved 
with zero initial conditions, P = P = 0att = 0. Define an auxiliary field Q by P = uj p Q. 
The Maxwell's equations are cast in the Schrodinger form: 

M = HV , (5) 

where the wave function \1/ and the Hamiltonian H are 



/ y/eE\ ( ice _1 / 2 Vx -i 



-IOJ. 



p 



V = B , H = -icV x e- 1 ' 2 . (6) 

\ Q / \ iu p -27 / 

The norm of the wave function, \\^\\ 2 — J dr^ty, is proportional to the total electromag- 
netic energy of the wave packet [23, 26]. When the attenuation is not present, 7 = 0, the 
Hamiltonian is Hermitian, and the norm (energy) is conserved. The quantities e and oj p are 
position dependent so that e — 1 everywhere outside the dielectric inclusions, and uj p = 
everywhere outside the metal part of the grating. 

It follows from Eq. (5) that ^(t + At) = exp(—iAtH)^/(t). In our simulations, the action 
of the infinitesimal evolution operator exp(-iAtH) on the wave function ^ is carried out by 
the algorithm described in [23]. The initial wave packet is Gaussian and propagates along 
the z axis perpendicular to the grating. Its spectrum is broad enough to cover the frequency 
range of interest. A change of variables is used to enhance the sampling efficiency in the 
vicinity of medium interfaces so that the boundary conditions are accurately reproduced 
by the Fourier-grid pseudospectral method [27, 28]. A typical size of the mesh corresponds 
to — l5D g < z < 15D g , and — 0.5D g < x < 0.5D g with 512 and 128 knots, respectively. 
The frequency resolved transmission and reflection coefficients are obtained via the time-to- 
frequency Fourier transform of the signal at some distance in front and behind the grating 
[29]. An absorbing layer is introduced at the grid boundaries in order to suppress artificial 
reflections of the wave packet [30] . 
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3 Results and discussion 



In Fig. 2a we show an interpolated image of the time evolution of the electric field E y along 
the z-axis passing through the center of the grating (see Fig. 1). The grating structure is 
characterized by e = 11.9, D g = 1.75 /im, a = 0.35/im and h = 1.4/xm. The red and blue 
colors correspond to positive and negative values of the field, respectively. The horizontal 
axis represents the ^-coordinate expressed in units of the grating period, D g . The grating 
extends from z = to z = h/D g . The vertical axis represents the propagation time measured 
in femto-seconds. The initial pulse impinging on the grating has a duration of approximately 
25/s. The instant when the pulse hits the grating followed by the main reflected signal is 
clearly visible in the figure. One also observes that a fraction of the electromagnetic energy is 
stored in the grating structure and leads to a long lasting radiation on both the transmission 
and reflection sides. This lasing effect extends to a pico-seconds time range, i.e., it is much 
longer than the duration of the initial pulse. It can be explained as due to the existence of 
trapped electromagnetic modes or resonances. It is the radiation of decaying trapped modes 
that comes with a phase opposite to the corresponding harmonic in the initially reflected 
pulse to the left from the grating structure and leads finally to the reduced reflection. The 
same lasing effect to the right from the grating structure is responsible for high transmission 
at the same frequency (see further discussion and results in Fig. 3). 

Figure 2b shows the dynamics of the electric field E y in the same setting, but the atten- 
uation 7 is set to zero. In this case, the trapped modes do not dissipate their energy into 
the metal. As a result, they live longer, which is clearly seen from comparison of the color 
intensity of the vertical strip in the middle of Figs. 2a and 2b (that represents the electric 
field of the trapped modes). 

The calculated transmission coefficient is presented in Fig. 3 as a function of the wave- 
length expressed in units of the grating period, D g . Each of the resonances observed in Fig. 
3 is associated with the corresponding trapped mode. The rough estimate given in Eq.(l) of 
their energies (frequencies) can be improved by taking into account the penetration of the 
field into the metal whose dielectric properties are described by Eq. (3). This yields 



(7i cm \ ^ 

h^) ' m = 1 ' 2 '-> ( ? ) 

where Q is the frequency of the lowest eigenmode in the stationary equation 

92E{X) +^e M (x,n)E(x) = 0. (8) 



dx 2 c 

Equation (8) is solved numerically under the condition that E(x) must decay exponentially 
outside the interval < x < a (i.e., in the metal). For a = 0.35 jim and h = 1.4/xm, 
this gives the following resonant wavelengths (expressed in the units of grating period): 
Ai = 1.493, A 2 = 1.352, A 3 = 1.185, A 4 = 1.031, A 5 = 0.899. The improved estimate of the 
resonant wavelengths agrees closely with the results obtained from numerical simulations for 
the exact problem only for the largest resonant wavelength (Ai = 1.493). The wavelengths 
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corresponding to the maximum of the transmission coefficient for higher modes are redshifted 
as compared to the estimated values. A similar result was reported by Takakura [7], but 
for the TM incident wave polarization. The red-shift can be explained by spreading of the 
trapped modes into the vacuum due to a finite thickness of the grating (see also Fig. 4 
and its discussion below), while the Fabry-Perrot modes of the electromagnetic field, used in 
our rough estimate, satisfy the zero boundary condition at the dielectric-vacuum interface. 
Clearly, an increase of the spatial volume occupied by a standing wave implies increasing its 
wavelength and, hence, lowering its frequency. 

Because of dissipative losses of energy in the Drude metal, the transmittance does not 
reach 100% and is, in fact relatively small. While for a lossless medium the sum of the 
reflection and transmission coefficients must be one as follows from the electromagnetic 
energy conservation, this is not the case when attenuation is present (blue curve in Fig. 
3). The maximal loss of energy corresponds to resonant wavelengths. This can be easily 
understood because the trapped modes remain in contact with the metal much longer than 
the main pulse (cf. Figs. 2a and 2b), and, therefore, can dissipate more energy through 
exciting surface electrical currents in the metal. We further illustrate this point by computing 
the transmission coefficient in the same system but without attenuation (7 = 0). The result 
is shown by the dashed black curve, which reaches 1 at resonant wavelengths. Observe the 
deviation of the transmission coefficient from 1 for the narrowest resonance at A ~ 1.5D g 
even in the absence of absorption. This resonance possesses an extremely long lifetime, so 
much so that we had to stop the calculation before it had decayed completely; that is, the 
total energy trapped into this mode was not completely radiated out and, hence, was not 
fully accounted for. For the resonance at A < D g the transmission coefficient does not reach 
1 because we study only the zero diffraction order scattering channel. 

Wave functions (field configurations) of the trapped modes can be extracted from the 
time-dependent wave packet by the time-to-frequency Fourier transform: 



The wave function of a particular trapped mode is obtained by setting 00 to the frequency at 
which the transmission coefficient attains the corresponding maximum. In order to improve 
the contrast, the wave functions have been extracted for the case with no attenuation. The 
results are presented in Fig. 4 for the grating with parameters e = 11.9, D g = 1.75 /im, 
a = 0.35/xm, and h = lA/im. There are five colored strips in Fig. 4. Each colored strip 
represents the electric field of the trapped mode associated with corresponding maxima of 
the transmission coefficient (see Fig. 3). The data (from the top to the bottom) starts with 
the largest wavelength resonance at A ~ 1.5D g and ends with the resonance at A ~ lD g . 
The red and blue colors represent, respectively, negative and positive values of the electric 
field amplitude. Each colored strip covers the coordinate range: —0.2Qfim < x < 0.26/zm 
(along the vertical axis) and the range for z (horizontal axis) is specified in the figure in units 
of D g . The fields localized inside the grating and fields radiated into the vacuum (the lasing 
effect) are clearly visible in the figure. The trapped modes localized inside the dielectric 
part of the grating exhibit a nearly Fabry-Perrot pattern with respect to the quantization 



00 




(9) 
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in the z-direction. Observe a slight spreading of the field into the vacuum regions, z < 
and z > h, which explains the redshift of the resonant wavelengths as compared to the pure 
Fabry- Perrot prediction give by Eq. (7). 

The structure of the field in the present case is such that one can regard the grating open- 
ings as an ensemble of independent emitters. They are coherently excited by the incident 
pulse, and their coherent emission builds up the radiation field associated with the resonantly 
enhanced transmission (reflection) properties of the grating. This is in contrast with the pre- 
viously reported TM results, where the excitation of plasmons leads to the coupling between 
effective emitters associated with the grating openings. This point is further illustrated in 
Fig. 5, where we show the transmission coefficient calculated for the grating structure with 
parameters e = 11.9, a = 0.3 fim, h = 0.8/im, and different periods D g = D = 1.75/im, 
D g = Do/1-5, and D g = D /A. The change of the grating period does not affect positions of 
the peaks in the transmission coefficient, pointing at the independence of the trapped field 
associated with different openings. On the contrary, for the TM polarized light, the resonant 
wavelengths are "pinned" to the grating period [5, 6, 11] - this fact being a reason for the 
ongoing discussion on the role of surface plasmons in the TM radiation transmission. The 
overall increase of the transmission coefficient when the grating period is reduced is due to 
the increase of the density of emitters (openings), while the time scale of the lasing effect 
remains the same because it is set by the attenuation of the metal and by the coupling of each 
individual grating region to the vacuum. In agreement with Eq. (1), a direct comparison of 
the result for D g = D with the results presented in Fig. 3 shows that the reduction of h, 
and primarily of a, leads to the blue shift of the whole resonance series. Similar dependence 
of the resonance wavelengths on the grating thickness h has been reported as well for the 
TM polarization [5, 6]. However, some caution is needed when comparing the TE and TM 
results in view of their different boundary conditions. 



4 The significance of trapped modes in resonant scat- 
tering 

The example of the grating structure considered above suggests that the knowledge of long- 
lived trapped modes is crucial for a custom design of nanostructured materials with enhanced 
transmission (reflection) properties in designated narrow intervals of wavelengths. Here we 
offer a rather general approach which establishes a direct relation between transmission 
(reflection) properties of nanostructured materials and the existence of trapped modes. We 
assert that once the Maxwell equations have been reformulated in the Schrodinger form, the 
significance of trapped modes for the light transmission can immediately be understood from 
the basic principles of quantum resonant scattering theory [35]. 

Recall that the approach relies on a representation of the total wave function of the system 
as a sum of non-resonant and resonant contributions [31, 32]. Consider the case without 
losses so that the corresponding Hamiltonian is Hermitian. Let H be the total Hamiltonian 
of a nanostructured material, supporting resonance(s), and H be a Hamiltonian responsible 
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for a non-resonant scattering. In the present case H and H are Hamiltonians of the metal 
grating with and without dielectric insertions, respectively. For the TM polarization the 
natural choice will be to set H to be a Hamiltonian of the simple metal slab with no 
gratings. The Lippmann-Schwinger formalism [33, 34] is applied to describe the scattering 
of a plane wave on a scatterer that has resonant excitations. Now we will show that the 
existence of long-lived quasi-stationary states in the symmetric [z — > —z transformation) 
grating with dielectric insertions implies that there exists a frequency at which the grating 
becomes transparent. If iv is the frequency of the incoming wave, then a solution of the 
Schrodinger equation (H — u)^/ = can be written in the form 

= + G + (u) (H - Ho) = + ^ + , (10) 

where G + (lu) = [uj—H+ir]]^ 1 , 77 — > 0+ is the Green function, \l>o is a solution of (H — iv)^f = 
0, and satisfies radiation (outgoing wave) boundary conditions and describes the scattered 
wave due to dielectric insertions. If the incoming wave is polarized along the grating, then 
the wave functions \l/o and contain only one component of the electric field which is 
denoted Eq and E + , respectively. Let the frequency iv be in the range in which the structure 
described by H Hamiltonian is a total reflector, then 

E -> A e ikz + A* e~ lkz , z -> -00 (11) 

and Eq vanishes as z — > 00 assuming that the grating is centered at z — 0. For the TM 
polarization, the magnetic field should be considered along similar lines. 

Quasi-stationary states, or resonances, that exist in the grating with dielectric insertions 
correspond to eigenvectors of H with outgoing wave boundary conditions and, therefore, 
they are associated with poles of the Green function G + {uo). Note that due to the complex 
boundary conditions H would have complex eigenvalues Uq — iT /2 with negative imaginary 
parts (r > 0) which specify resonance widths. Thus, in the vicinity of a pole, the frequency 
dependence of the Green function can be approximated by 

G + {uj) 1 — . (12) 

Let iv be near a resonant frequency (within the resonance width), while the resonant fre- 
quency c^o is assumed to remain in the range of total reflectivity of the pure metallic grating. 
In this case, from the symmetry of the Hamiltonian H under the parity transformation, 
z — > — z, it follows that 

E + -> [s\gn{z)] p A + e ±lkz , z -> ±00 , (13) 

where the parity factor of A + corresponds to either symmetric (p = 0) or anti-symmetric 
(p = 1) eigenfunction of H. The lowest frequency resonance corresponds to the symmetric 
solution. The amplitude A + has to be found from the energy flux conservation. The incident 
flux is I ^4o 1 2 - The outgoing flux is |Ao±A + | 2 + |A + | 2 . Let 0o and <p + be the phases of Aq and 
A + , respectively. Then from the flux conservation we infer that \A + \ 2 = \A \ 2 cos 2 (0 o — 4> + )- 
From Eq. (12) it follows that + , as a function of the frequency tv, rapidly changes by ir over 
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a small interval containing the resonant frequency (the eigenvalue of H), while the phase <J)q 
describing the non- resonant scattering is nearly constant, or changes slowly. Therefore by 
continuity of <j) + in the vicinity of the resonant frequency there exists a frequency at which 
\A \ 2 = \A + \ 2 , that is, the incoming flux coincides with the transmitted flux and the grating 
becomes transparent. 

Thus, in the absence of attenuation, the existence of trapped electromagnetic mode(s) or 
resonances necessarily leads to 100% transmission close to the wavelengths of the trapped 
modes in symmetric nanostructured materials which otherwise (in the absence of such modes) 
are not transparent. The same approach can be used to analyze a possible 100% reflection in 
dielectric gratings, which otherwise (in the absence of trapped modes) are nearly transparent 
in the zero diffraction order. It also explains possible Fano profiles in the transmission (re- 
flection) coefficient in the cases when non-resonant scattering described by the Hamiltonian 
H leads to both reflection and transmission. 

In a generic case, in order for 100% transmission to be possible, the parity symmetry of 
the Hamiltonian (system) is not required, but the weaker condition (13) on the asymptotic 
behavior of is indeed necessary. The latter readily follows from the flux conservation. A 
sufficient condition for 100% transmission is the absence of scattering channels with different 
quantum numbers at the resonance frequency, e.g., higher order diffraction and/or coupling 
to scattering states with different polarization. All these effects lead to breaking of the 
flux conservation in a selected scattering channel. If the coupling to other channels is not 
significant, an enhanced (not 100%) transmission can still be observed in the selected channel. 
It can be quantified by conventional means of quantum scattering theory applied to the 
effective electromagnetic Hamiltonian of the system in question. In practice, gratings often 
have the imperfections which induce a coupling between TE and TM polarization. Thus, 
in the TE resonant scattering channel 100% transmission will be lost due to leaking of the 
energy flux into higher order diffraction and/or TM polarization channels induced by the 
imperfections. One should not, however, expect a resonant transmission in the TM channel 
facilitated by the imperfections since the resonant (enhanced) transmission is essentially 
due to the constructive interference which, in turn, occurs thanks to the periodicity of the 
structure, while imperfections are usually randomly distributed. In fact, in a lossless grating, 
a deviation of the transmission coefficient from one at the resonant frequency can be used 
as a measure of the grating quality. 

Finally, it is worth mentioning that a loss of the electromagnetic energy in dispersive 
materials prevents 100% transmission (reflection) to occur, again because of breaking of the 
flux conservation even in a single scattering channel available. This is clearly seen from our 
numerical results with and without attenuation (Fig. 3). Our formalism offers a possibility 
to quantify such effects by studying the unitarity violation of the scattering matrix caused by 
the skew-Hermitian part of the effective Hamiltonian (see, e.g., [36]). Note that in quantum 
systems a leak of the probability density into scattering channels weakly coupled to the 
one of interest is often modelled by an effective non-Hermitian Hamiltonian for the main 
(selected) scattering channel(s) only. 
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5 Conclusions 



We have elucidated the role of trapped modes in the extraordinary light transmission in 
nanostructured materials by reformulating the Maxwell equations for passive linear media 
in the form of the time dependent Schrodinger equation and applying to the latter the basic 
principles of quantum resonant scattering theory, in particular, the Lippmann-Schwinger 
formalism. Trapped electromagnetic modes in nanostructured materials play the same role 
as resonances in quantum scattering. This offers well developed quantum mechanical tech- 
niques to study resonant light transmission and reflection properties of gratings and other 
nanostructured materials. We have illustrated this approach by a detailed numerical study 
of a metal-dielectric grating. In particular, for the TE polarization of the normal incident 
radiation, the grating, while being a total reflector in the zero diffraction order when no 
dielectric fillings are present, has been shown to become transparent for certain (resonant) 
wavelengths when the fillings are present. In accord with the quantum resonant scattering 
theory, stationary states have been observed in the latter case and none in the former. 
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Figure captions 

Fig. I. Schematic representation of the studied system. The radiation is incident 
along the normal to the slab (the z-axis), consisting of a grating structure with alternating 
regions of metal and dielectric along the x-direction. The dark grey and shaded regions 
correspond, respectively, to the metal and dielectric parts of the structure. The system is 
translationally invariant along the y-direction. 

Fig. 2a. (Color online) Interpolated image of the time evolution of the electric field 
E y along the z-axis passing through the center of a dielectric region of the grating (see 
Fig. 1). The grating structure is characterized by e = 11.9, D g = 1.75 jum, a = 0.35/zm, 
and h = 1.4/xm. The red and blue colors correspond, respectively, to positive and negative 
values of the field with the color intensity related to the field magnitude. The horizontal axis 
represents the z-coordinate expressed in units of the grating period, D g . The grating extends 
from z = to z = h/D g . The vertical axis represents the propagation time measured in 
femto-seconds. The metal is described by the dielectric function of Eq. 3. 

Fig. 2b. (Color online) The same as Fig. 2a, but with no damping inside the metal 
(7 = 0). 

Fig. 3. (Color online) Zero-order transmission coefficient as a function of the wave- 
length of the incident radiation measured in units of the period D g . The calculation is 
carried out for the grating structure characterized by e = 11.9, D g = 1.75 /im, a = 0.35/xm, 
and h = 1.4/im. The solid red and dashed black curves correspond to calculations done 
with (7 = O.leV) and without (7 = 0) damping in the metal, respectively. The sum of 
the reflection and transmission coefficients for 7 = O.leV is shown as the blue curve. Its 
deviation from 1 represents the loss of electromagnetic energy because of the absorption in 
the metal. 

Fig. 4. (Color online) Electric field of the trapped modes inside the grating as function 
of x and z coordinates. The results are presented for the grating with parameters e = 11.9, 
D g = 1.75 /im, a = 0.35/xm, and h = lA/im. Each colored strip represents the electric field 
of the trapped mode associated with corresponding maxima of the transmission coefficient 
(see Fig. 3). The data (from top to bottom) starts with the largest wavelength resonance 
at A ~ 1.5D g and ends with the resonance at A ~ lD g . The red and blue colors correspond, 
respectively, to positive and negative values of the field. The data has been normalized to 1 
at maximum so that color scale covers the range [—1, +1]. The rr-range for each strip (along 
the vertical axis) corresponds to —0.26fim < x < 0.26fim and the range for z (horizontal 
axis) is specified in the figure in units of D g . 

Fig. 5. (Color online) Zero-order transmission coefficient as a function of the wave- 
length of the incident radiation measured in units of the period D = 1.75/im. The 
computed data are given for a grating structure characterized by e = 11.9, a = 0.3/im, 
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h = 0.8/xm, and different periods D g . The red curve shows the transmission coefficient for 
D g = D = 1.75/im, the blue curve for D 9 = D /1.5, and the black curve for D g = D /4. 
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